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Abstract 



In this paper we study a one-dimensional space-discrete transport equation subject to additive Levy 
forcing. The explicit form of the solutions allows their analytic study. In particular we discuss the invari- 
ance of the covariance structure of the stationary distribution for Levy perturbations with finite second 
moment. The situation of more general Levy perturbations lacking the second moment is considered as 
well. We moreover show that some of the properties of the solutions are pertinent to a discrete system 
and are not reproduced by its continuous analogue. 
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1 Introduction 

In what follows we discuss the system of coupled linear stochastic differential equations related to the one 
considered in Rcf . |MSVE07a] , which after time-integration reads 

'<(i;s) = 0, 

^ a'^{t;s) = ai{s) + / {-a^iu; s) - iya'^{u; s)) du + Ls{t), /-.x 



with 1/ > and for — oo<s<t<oo. The discussion in [MSVEOTa] was motivated by a description of 
hydrodynamic turbulence within a framework of simplified shell models, see |MSVE07b] . The system ([T]) is 
a space-discrete transport equation with a Dirichlet boundary condition and a singular random perturbation 
at the position n — 1. The interpretation of the equations corresponds to the situation when the energy is 
pumped into the system at shell n = 1 via the effective noise (or signal) term La{t), is transported between 
the shells and dissipated due to the effective viscosity v >Q. 

Let L = {L{t))t>o be a one dimensional Levy process, i.e. a stochastically continuous process with 
stationary independent increments starting at the origin and having right-continuous sample paths with left 
limits, and 




'ri+l 



(u] s) — va^-^{u; s)) du, n>2, t>s, 



Ls{t) ■.^L{t-s)-L{s), 



GO < S <t, 



(2) 



be a time-shifted Levy process defined on the half-line [s,+oo). It is well known that a Levy process L is 
completely characterised by the Fourier transform of its marginal distributions which has the following form 
due to the Levy-Hinchin formula: 

*(A) - -cj'^ + zAm - / (e'^'^ - 1 - ^A-^) p{dy). 

The function ^'(•) is called the cumulant of L. The first summand of '5 corresponds to a Brownian motion 
with variance cr, cr > 0, the second summand determines the linear drift ^f, £ R, and the third summand 
corresponds to a pure jump Levy processes whose jump intensity and sizes are governed by the jump measure 
p satisfying the integrability condition /ij\^|gj(2/^ A 1) p{dy) < oo. 

In particular, if p = 0, the Levy process L is just a Brownian motion with drift. The case of the standard 
Brownian motion without drift, corresponding to the cumulant 5'(A) = — A^/2, is exactly the one discussed 
in Ref. [MSVEOTa] . Further, if p(R\{0}) = c g (0, oo), and tr = /x = 0, then L is a compound Poisson process 
with intensity c and jump sizes distributed with the probability law p{-)/c. Finally, the case a = p = and 
p{dy) = c{a)\y\~^~°'dy, a S (0,2), corresponds to Levy flights with ^'(A) = — |A|". We address the reader 
to the books [Sat99[ |App04| for more information on Levy processes. 

The goal of this paper consists in studying the long time behaviour of solutions of the system ([T|) . Our 
approach here differs from the one of |MSVE07a] and is based on the explicit solution of ([T]) , which allows a 
generalisation of the results of |MSVE07a) . We moreover discuss some peculiar properties of the continuous 
limit of ([1]) described by a partial differential equation. The behaviour of the solutions of this continuous 
equation turns out to be different from the one of its discrete analogue, Eq.(l). 

2 Explicit solution 

Let us first derive the explicit form of the solution of ([T|). 
Proposition 2.1 The solution of ([T]) is given by 
a^(i;s) = 0, 

OO 

a^i; s) = e-''(*-^) °™(«) [j\n-m\m - s)) + (-1)—! J„+„(2(t - s)) 

m—1 

rt (4) 

-t- / Hi;{t-r)dLs{r), n>l,t>s, 



J„(-) being the Bessel functions of the first kind. 

Proof: For brevity we set s = 0, aj^ a„(0) and aj'j(t) := a^{t\Q). Since the system ([1]) is linear it is 
enough to solve the homogeneous deterministic system for initial conditions aj^ = 1, m > 1, a^^ = 0, n ^ to, 
and the non-homogeneous stochastic system with initial conditions aj^ = 0, n > 0. Thus the solution of the 
homogeneous system is given by 

<(t) = e-'^*(J|„_„|(2i) + (-l)'"-ij„+™(2i)), 
<nit) = e--*(Jo(2t) + (-1)™-V2™(2t)). 

The proof is straightforward with help of the differentiation formulae (here and below, see [AS84| Chapter 

9]) 

24(t) = J„_i(t)- J„+i(t), n>l, 

4it) = -Mt). 

The initial condition is satisfied due to the relations Jo(0) = 1, J„(0) = 0, n > 1. 
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To show that H!^{t — r)dL{r) solves the non-homogeneous equations we have to be sure that the 
stochastic integral is well-defined for all t > 0. First we note that the asymptotic expansion 



J„(2r) « — , 



0, n > 1, 



(7) 



implies that the integrand is a smooth and bounded function on r S [0, oo) for all v >Q. Consequently the 
stochastic integral is well-defined and we can easily calculate the characteristic function 



Eexp(iA^ Hi^it - r) dL{r)^ =exp(^ ^'(AiJ,';(r)) rfr) , Ael 
For the next step of the proof we take use of the relations ^ to show that for n > 1 



dr 



Kir) 



2r2 

J„-i(2r) + J„_i(2r) 

Ys ' +' 



Jn-i{2r) - J„+i(2r) -nj„(2r) ^_^^ 



2r 

Jn-i{2r) - Jn+i{2r) - uJ^{2r) 
2r 



(8) 



(9) 



Thus wc finish the proof with help of the intergation by parts. Indeed for n > 1 we have 

ft ^ fU 



Hi:it^r)dL{r)^Hi:{0)L{t) + 



m{t-r)dL{r) 



du 



H::iO)L{t) + / - - ^^<.H) du. 



(10) 



Taking into account the asymptotics ([7]) we note that H^{0) = 1 for n = 1 and H'^^{Q) = for n> 2, what 
finishes the proof. ■ 



3 Stationary distribution 

In this section we study the convergence of a'^{t, s) to the stationary distribution. 

It is clear that due to the linearity of the system, the convergence to the stationary regime depends on 
the asymptotic properties of the deterministic solution corresponding to L = 0. In [MSVEOTa] . the authors 
showed that for (a'^)n>i S > 0, the deterministic solution converges weakly to zero. For bounded initial 
values of a„, they gave an example of no convergence, namely the solution 

aL 1 W = 1, 

n (11) 
aLW = 0, n>l, 

which is a fixed point of the homogeneous system. With help of the explicit formula for the solution, we 
can find another bounded initial condition, for which the corresponding solution weakly converges to zero. 
Indeed, for a2n-i = 0, 02,1 = 1, n > 1 we have: 

n-1 

n-1 (12) 
OanW = M'^t) + 2 J2 ^2j(2t) + J2n(2t). 

The convergence follows from the asymptotic expansion 

Jn{2t)^ ^ — cos[2t- — - -j, t^oo. (13) 
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The convergence is however not uniform over n > 1 since for any t > 

i.2t 



2Vj2j„i(2i)^ / Jo{r)dr 



(14) 



and 

n-l 



Jo(2t) + 2^ J2j(2i)^l, oo. (15) 

To study the influence of the stochastic forcing, let us set aj^{s) := for n > 0, so that the deterministic 
part of the solution disappears. To obtain the invariant distribution, we consider the pull-back limit 

pO poo 

al= hm / H::{-r)dL,{r)^ K{r)dL{r). (16) 

s^-^Js Jo 

It is instructive to determine the Fourier transform of a'^ which can be found as 

F,^(A)= lim Ee*^<("'") =Eexpf / '^{XHi^ir)) dr) , A S R. (17) 

First of all, wc have to study the existence of the latter limit. Clearly, is a well-defined random variable 
if and only if for any finite A 



1'(Aff;;(r))dr < oo. (18) 
In view of the asymptotics this is equivalent to the convergence of the integral 

'^/e-cos(2r)x ^^^^ (19) 



which in turn depends on the asymptotics of ^'(r) as r ^ 0. We also notice that the condition (fTO|) does 
not depend on n, thus all random variables a"^ are either well-defined or not well-defined simultaneously. In 
particular, for any n, m > 1 we can calculate the mutual Fourier transform 

K,ni>^n, A„) := Ee'(^"<^+^"<-) - exp ( - ^{XnK{r) + X,nK,{r)) rfr) . (20) 

Clearly, similar formulae hold for any finite number of summands. 

In the next section we study the invariant distributions a'^ in more detail by discussing several special 
situations corresponding to difi'crcnt forms of the cumulant 'I'. 



4 Particular cases 
4.1 Linear perturbation 

We start with the simplest deterministic perturbation L{t) = ^it linear in time (corresponding to a constant 
forcing) with the cumulant 5'(A) = iA/i. The Fourier transform of the stationary solution is found explicitly 
as 



Ee*-^"" = exp ( ^ «'(Aij;;(r)) dr 

= exp{iXtiJ K{r)dr) 

( p J„(2r) 
= exp I zA/i / n e dr 



U2r) ■ 



exp (iA^77,( ^ + \/ 1 + ^ 
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This means that the stationary solution is a constant 



< = /^"(^ + yi + ^) " (22) 

and as — !• wo obtain that 

a'^ K ^in(l - ^-^ly + . . ^ ^ ^j,n. (23) 
On the other hand as n ^ oo, aj:^ ^ for i/ > 0, whereas a° — > oo. 

4.2 Gaussian perturbation 

If the forcing is Gaussian, i.e. ^'(A) ~ — it is clearly seen that each a^^ is also a Gaussian random 
variable with the characteristic function 



I]e'^< = exp ( ^ *(AiJ,';(r)) dr) 

.(-A^^/V;;,.,)^*) 



exp I 



^,1(2r)_ 



exp -A"V / n':^^^e~"''dr 



n 



A„=A„=0 



2/o..xi-m-n (m + n-2)! 



(m- l)!(n- 1)! 

TO + n— 1 m + ri m + n + 1 m + n , , ,4 



/TO + n— i m + ri m + n + i m + n 
X aF^\^ , — - — , , — h l;m + l,n + l,m + n + 1; 



,•2 



^nn 92 / 771— n\ 

Ea^a'^ = ^ - cos J 



1 1 



(24) 



/ A2cr2r(n- i) , o„ / 1 1 4 \N 

= exp ( - ^^^'^'-'-sF^ (n - n, . + -; n + 1, 2. + 1; --) ) . 
In the limit as ^ we recover 

Ee'-^ - exp f - ^(^^^ -n. + ...))^ Ee-< = exp f - ^ f 1 - JL)") (25) 



2 V7r(4n2-1) 'V tt V in^ J 

for all 71 > 1. In particular, 

E«)2 = n^oo, (26) 

as it was discovered in |MSVE07a] . Taking into account Ea. ipO)) . we determine the covariances, 

a2 roc 



(27) 



- (to + 7i)2 — 1 (rn — nY — 1 
where the latter formula consides with the result from |MSVE07a] . 

4.3 Levy perturbations with finite second moment 

Further, we note that the covariance structure of the stationary distribution is the same for all Levy forcings 
with the finite second moment. Indeed, assume that EL(t)2 < 00 which is equivalent to the condition that 
^'(•) G C2(]R). Differentiating the Fourier transform (|20p at zero, we determine the mean value and the 
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covariances of the stationary solutions aj:^: 

d 



E< 



EU V 



dX 



A=0 



*vl''(0) / F:(r)dr, 



F:{X)\^^^ = -^"{0)J^ {K{r)fdr-{^'{Q)j^ K{r) dr) 



92 



n 



(28) 



A„=A„=0 



= -*"(0) / Hi;^{rmr) dr - (vI/'(0))2 / H'^{r) dr / H'^{r) dr. 
Jo Jo Jo 

The closed form for the above integrals was already given in (HT|). (pij) and ([TT)) . Thus our analysis shows 
that the covariance structure of the stationary laws does not depend on Gaussianity itself, and is the same 
for all Levy forcings with equal first and second moments. 

4.4 Levy flights 

Finally, for the Levy flights forcing with ^'(A) = — (t"|A|", a G (0,2), a > 0, the stationary probability 

distribution exists for all n > 1, a £ (0, 2), and > and for all a S (f 7 2) and u = 0, and is also a-stable 
with the characteristic function 



Ee*^°" = exp I 



dr). 



(29) 



In the zero viscosity case 1^ = 0, we observe a strange transition at the critical value a = 2/3 of the stability 
index. Formally it comes from the integrability condition of the Bessel functions. Indeed, in view of the 
asymptotics ([T3| the limits 



f*\Jj2rT 
lim / dr, iim 



icos(or 

„3a/2 



dr and lim 



dr 



(30) 



either exist or not exist simultaneously, and the latter limit is finite if only if a > 2/3. Since the stability 
index determines the weight of the tails of the Levy flights, P(|i(t)| > u) = ©(it^"), u 00, we can 
conclude that in the model under consideration, too heavy tails and big jumps of the random forcing lead to 
anomalous dispersion of the energy in space and do not allow to build up a stationary probability distribution. 
On the other hand as it will be shown in the next section, the critical value a = 2/3 has it origin not in the 
nature of the random forcing but rather in the properties of the space-discrete deterministic equation. 



5 Continuous limit 

In this section we discuss an interesting effect related to the continuous limit of Eq. U]) . For any real function 
/ we introduce a discrete symmetric gradient at x as 



V,J(x) 



fix + h) - fix ~ h) 
2h 



h>0. 



(31) 



Denote also a^it, n; s) := a'^it; s). With this notation under condition h — 1, the system ([T]) takes the form 
'a''{t,0;s) = 0, 

a'^ {t, h] s) = a{h; s) + J i—2\/iia'^iu,h;s) — i^ha'^{u,h;s))du+j^Lgith), ^^2^ 
a'^{t,nh;s) = a{nh;s) + / i—2'S/ha'^iu,nh;s) — i'ha'^iu,nh;s))du, n>2, 



V/i denoting the discrete gradient w.r.t. the second variable of a'^ . We omit the discussion on the limiting 
behaviour of the deterministic solution assuming that the initial values of a" are identically zero. According 
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to (|16p . we find the stationary solution of the perturbed system as 

J-s -r " hj„ 



= hm / ° -^"^ ^'^K '^'^- dL,{hr) - ^ ^^^^e^'^''' dL{hr), (33) 

-oo ft ' - t, I „ 



provided the latter stochastic integral exists. To study the existence, we use the change of variables formula 

f{r)dL{hr)l f{i)dL{r) (34) 
Jo 

for h > and a continuous function / to obtain that 

u / \ d n f ^ri(2r) _,^f^r IT /u \ d, n f Jn{ — ) jtI \ 

ahin) = T / e dL{hr) = - / — dL[r). (35) 



It is clearly seen, that the stochastic integral on the right-hand side of ((35)) exists for all n > 1 and /i > if 
and only if the integrability condition (fT9|) holds (recall constraints on a in the Levy flights case). Further, 
let n and h be such that nh ^ x > Q. Then with help of (|35p we can pass to the limit 



aY{x) ^ r "-4" r ^AI ■dL(r) ^ -er^^ l(-) =: a''{x). (36) 

Jo h r Jo r x \2J 

To derive the latter limit we also used that (fiJ^{^{x + l-)))p,yo, is an approximating sequence for the Dirac 
(5-function (see |Lam69b] ). i.e. that for any appropriate test function / we have 

/•oo 

lim / fiJ^Ui{x + l))f{x)dx^ f{0). (37) 

This approximating system which appears here naturally as a generic property of the model, also proved 
to be useful in treating relativistic plasma dispersion in the limit of zero magnetic field, see }Lam69aj . We 
emphasise that the limit pop is pointwise for a; > and /i — *■ 0, and our argument has sense provided the 
condition ([TO]) holds and the random variables a'j[{x) are well-defined. 

As a partial differential equation analogue of (j32[) . let us consider the linear transport equation 



b^'^it, X- s) + 2Vb''^^{t, x; s) + vb^^^t, x; s) = L,{t)5{x - e), t > s, a: > 0, 
b^'^is.x-.s) = if{x), x>Q, </5(0)=0, (38) 
6'''^(f,0;s) = 0, t>s, 

where > 0, e > 0, (5 denotes the Dirac function, and tp is the initial condition. Such a setting allows to 
capture the random perturbation as a non-homogeneity in the linear transport equation. 
Equation (j38|) can be easily solved explicitly. Its general solution is given by 

b''-^{t,x;s) = <^{x-2t)+ I e^''^'-'''>d{2r-2t + x-e)dLs{r) (39) 



for an arbitrary <&(•). The first boundary condition yields: 

<^{u) = ip{u + 2s), u > -2s. (40) 

The second boundary condition yields: 



$(-2i) = - / e-''^'~'''>6{2r - 2t - e) dLs{r), t > s, 
$(u) = - / ' e''(t+'')(5(2r + u-£)rfLs(r), u < -2s. 



Consequently, 

b'''''{t,x;s) = ipix-2t + 2s)+ I e'"'^*-''"' 5{2r ~ 2t + x - e) dLs{r), s<t<s + x/2, (42) 
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and 



b''^t,x;s) = j e-''^'-''H{2r - 2t + x) dL{r - s) - j e-"^'"-'^^ S{2r - 2t + x - e) dLsir), t>s + x/2, 

(43) 

In particular, for t = 0, fixed a; > and s < — x/2 we have 

f \^'''-^+^^S(2r + x-e)dLjr) 

(44) 



b^^'it^x^s)^ f e''''S{2r + x-e)dLs{r) - I ' e''^''+^^S{2r + x - e) dLs{r) 



= e-'''^L{-^) e-'^tL(-f ) =: f'ix). 
This yields the stationary solution for the limiting equation for v > 0: 

6-(x) =e-''5L(|), a;>0. (45) 

It is important to notice that in this case we have no constraints on the forcing L. 

As we see, the stationary distributions of a'^{x) and b''{x) are different due to the pre- factor ^. However, 
for i> > the influence of this pre-exponential term is not crucial at least for large values of x. Another 
difference concerns the fact that the asymptotic solution a"^ can be obtained only for a class of forcings 
satisfying condition (flQl) . whereas the solution 6"^ is defined for all Levy drivers L. The reason for these 
differences can be the inaccuracy of the approximation of the discrete gradient V/i by gradient operator 
which neglects higher order derivatives. Due to the presence of higher derivatives in V?i, oscillating solutions 
of the deterministic system arise in a discrete system. These lead to complex interference between the 
modes in the discrete model. This interference makes a'^ to be typically smaller in amplitude than the 
solution of the continuous model where no such complex interferences arise. Moreover, the interference 
between many modes in the case of heavy-tailed forcings leads to the divergences and implies non-existence 
of the corresponding stationary probability distributions. Technically, the reason for this can be explained 
as follows. The random perturbation comes in as a stochastic integral of an oscillating kernel J„(2r)/r which 
under certain conditions can be approximated by the (5-function. However the asymptotics of this kernel 
at infinity imposes restrictions on the space of admissible test functions, as it happens in the case of zero 
viscosity and Levy perturbations with heavy tails and, i.e. when v = and P(|i(t)| > u) = 0{u~°'), u — > oo, 
a G (0,2/3]. 



6 Conclusions 

In this paper we studied a space-discrete transport equation subject to additive Levy forcing. We derived an 
explicit solution of this equation in terms of the Bessel functions of the first kind, and showed that the random 
forcing comes into solution via a stochastic integral of a certain oscillating kernel. We generalised the results 
by [MSVEOTa] on the covariance structure of the stationary distribution from the purely Gaussian to the 
general Levy case, and in particular established that the covariances are identical for all Levy forcings with 
equal first and second moments. In case of zero viscosity and non-Gaussian Levy flights forcing, we found 
that the stationary probability distribution exists only for stability indices a £ (2/3,2). The explanation 
to this phenomenon is given by comparison of a space-time limit of the space-discrete transport equation 
with its partial differential equation analogue. An interesting approximating sequence for a Dirac (5-function 
appears as a natural part of the space-discrete system. In particular, this approximating sequence allows 
only for test functions which have tails lighter than u~^l^ . The discrepancy between the space-discrete 
and the PDE models comes from taking into account higher derivatives while considering a space-discrete 
gradient operator. 
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